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ABSTRACT 

Theory has been developed for predicting the loss of 
response of a QCM to a liquid deposit due to viscous effects 
in the deposit. The Iojs of response is expressed by a 
response factor, equal to the response of the QCM to a liquid 
film divided by its response to a solid film of the same mass 
per unit area. The theory assumes a droplet-type deposit 
morphology, and considers the influence of droplet distribution 
parameters. Experiments have been conducted to examine the 
validity of the theory, using DC 704 silicone oil as the subject 
deposit material. Experiments were made in two series - one 
with constant deposit mass and variable temperature, the other 
with variable deposit mass and constant temperature. Satis- 
factory agreement with the theory was found. Interpretation of 
the data using the theory has enabled information on droplet 
area coverage and number density to be deduced. 

INTRODUCTION 

The quartz crystal microbalance (QCM) has been used in 
the aerospace industry for several years to monitor the build- 
up of condensed contaminant deposits on spacecraft surfaces in 
Doth simulation and flight situations. Also, the QCM is being 
used increasingly in basic research on contamination- related 
phenomena and its accuracy and behavior in this application 
have consequently been subject to critical analysis. One of the 
areas of concerr is the response of the QCM to liquid deposits, 
as opposed to solid deposits. Because the dynamic coupling of 
a liquid deposit to the crystal surface is less than one hundred 
per cent, the sensitivity of the QCM, measured in frequency 
change per unit added mass per unit area, is reduced. This is 
a significant problem in the aerospace contamination measure- 
ment application because many of the species outgassed from 
common materials, e.g. RTV silicone adhesives and potting 
compounds, as well as common diffusion pump oils exist as 
very low vapor pressure liquids at near ambient temperatures. 
It is therefore desirable to develop a technique for predicting 
the response of the QCM to liquid deposits quantitatively. 

When a liquid deposit is oscillated in shear by the QCM 
measuring crystal only the liquid layer in immediate contact 
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with the crystal moves with the same velocity as the crystal 
surface. Liquid layers some distance from the surface can 
move only in response to a shear stress in the liquid, which 
requires a velocity gradient to be established. The presence 
of a velocity gradient means that the average deposit velocity 
during oscillation will always be less than that of the crystal 
surface. It has been shown elsewhere'^) that the frequency 
change of the QCM produced by an added mass is a function of 
its kinetic energy. Since a solid deposit moves throughout with 
the same velocity as the crystal surface, the change in fre- 
quency produced by an added liquid mass will always be less 
than that produced by an added solid mass of the same magnitude. 

A further reduction in QCM response is produced by the 
tendency of liquid deposits to occur as a distribution of droplets. 
The formation of a droplet from an equivalent mass uniformly 
distributed over the same base area will increase the mean 
distance from the deposit mass to the crystal surface. Since 
the velocity of the liquid deposit decreases with distance from 
the crystal surface the formation of droplets will decrease the 
average velocity of the deposit, and hence further reduce the 
magnitude of the induced frequency change. The response will 
also be affected by the droplet size and number per unit area. 

For a given deposited mass per unit area the response of a 
system of a large number of small size droplets will be greater 
than that of a few large size droplets, because of the greater 
average distance of the deposit from crystal surface in the 
latter case. Also, the measured mass per unit area will vary 
with the number of droplets per unit area for a given droplet 
size. 

The approach taken in this paper to estimate the magnitude 
of these effects is to determine the velocity profile of the oscil- 
lating liquid deposit, and then to derive the QCM response 
using the Rayleigh energy method for harmonically oscillating 
systems. The problem is addressed theoretically in three 
stages. First, the effect of viscosity is determined for a uni- 
form liquid deposit. Second, the effect of droplet geometry is 
derived. Third, the effect of droplet size and distribution is 
investigated. A limited experimental program was then con- 
ducted to assess the validity of the derived theory. 

QCM Response 

The sensitivity of a QCM, S, is defined as the relationship 
between the added solid mass per unit area, (Am/ A) , and the 
resulting change in frequency, Af o . Hence for a soli% mass 

Af g = S g x (Am/A) g (1) 

This paper is not concerned with the absolute sensitivity of 
the QCM to liquid deposits, but with the relative loss of 
sensitivity if the deposit is liquid rather than solid. This can 
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be expressed by a "response factor", F, defined as the fre- 
quency change, induced by a liquid deposit &f , divided by the 
frequency changed induced by a solid deposifof the same mass 
per unit area. Hence, for liquid films the relationship equiv- 
alent to (1) is 

Af^ - S^x (Am/ A)^ (2) 

and, for (Am/A) = (Am/A) 
s L 

F •S 1 /S. = tfj/if, (3) 

F can be derived using the same Rayleigh energy method for 
undamped harmonically oscillating systems as was used 
success fully (1 ) (2) to derive the absolute sensitivity, S. In 
this method the total kinetic energy, T, and potential energy, 

V, of the system are expressed in terms of the displacements, 
u., in each degree of freedom, as follows 

T = ia(u 2 + u 2 %) (4) 

V = ic (u 2 + u 2 u 2 n ) (5) 

The free ency of the system, f, is then given by 

f 2 = c/a (6) 

For small perturbations of the system, such as the addition of 
mass tc the QCM crystal, equation (6) can be rewritten 

(f + Af) 2 = (c + Ac)/(a + Aa) (7) 

The surface of a crystal oscillating in thickness shear is an 
antinode, so mass added there does not affect the potential 
energy of the system, making Ac zero. Ignoring squares of 
small quantities, equation (7) can be rewritten as 

Af/f = -Aa/a (8) 

Equation (8) can be written for solid and liquid deposits. Noting 
that f and a are independent of the phase of the deposit, the 
following expression is obtained for F from (3) and (8). 

F = Aa /Aa (9) 

l s 

The response factor is thus determined from the ratio of 
kinetic energy coefficients of the liquid and solid deposits. 
These coefficients are determined from an equation similar to 
(4). The oscillating crystal and deposit is a single degree of 
freedom system, so the velocity of the system can be repre- 
sented by the velocity at a single location. The logical location 
to select in this context is the crystal surface. In general, 
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oscillation velocity and deposit density can vary across the 
crystal and these variations can be represented by a distri- 
bution function. (2) The effect of non-uniform distribution 
functions is to introduce a multiplicative constant into the 
expression for T and hence a. Since the present analysis is 
icerned with comparing equivalent situations this constant 
uould appear in both denominator and numerator of the right- 
hand side of equation (9), making F independent of the distri- 
bution functions. In this paper the distribution functions are 
thus assumed equal to unity to minimize algebraic complexity 
of intermediate expressions in the following derivations. A 
general expression for the kinetic energy of a deposit,AT, can 
then be written in terms of a uniform crystal surface velocity, 

v sin 2nft 
o 

AT = ^Aa(v sin 2nft) 2 (10) 

o 

The problem is now to obtain expressions for the kinetic 
energies of the liquid and solid deposits. This is done by 
integrating the kinetic energy throughout the deposit, for which 
the velocity distribution in the deposit must be known. The 
solid has the same velocity throughout, equal to that of the 
oscillating crystal surface. The velocity in the liquid deposit 
will vary through the thickness and must be determined by 
viscous theory. 

VELOCITY DISTRIBUTION IN THE DEPOSIT 


Uniform Film Deposit 


Figure la shows the model for the oscillatory motion of a 
uniform liquid film. The fluid motion in the film is described 
by the Navier-Stokes equations. For the present case there 
are nc pressure gradients or body forces, and flow takes place 
only in the y-direction. The Navier-Stokes equations thus 
reduce to: 


o v _ . 5 2 v 

3t " p ax 2 


(ii) 


This equation is to be solved for the boundary conditions: 

x = 0, v = v sin 2nft (crystal surface 

velocity) 


x = x , ^ = 0 (zero shear stress at 
111 x outer boundary) 

Equation (11) and the above boundary conditions are whol 
analagous to the heat conduction problem of a parallel siued 
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infinite slab with harmonically varying temperature at one sur- 
face with the other surface insulated. Carslaw and Jaeger give 
the following solution to this problem^), translated here into 
momentum transfer terms. 


cosh k(x -x) (1-ri) 
rn 

cosh kx (l+i) 
m 


/cosh k(x -x) (l+i) 
m 

\cosh kx (l+i) 
m 



( 12 ) 

(13) 

(14) 


v/v q = B sin (2nft+0 ) (15) 

B is the amplitude of the oscillation at a distance x from the 
crystal surface while 0 is the phase angle by which the dis- 
placement at this location trails that -\t the crystal surface. 
Numerical values of B and ) are given in Figure 2 as functions 
of kx^ and x/x^. Equation (12) can be expanded as follows: 

v/ v = Bcos0sin2rrft + Bsin0cos2nft (16) 

o 

The velocity of the liquid thus consists of a component Bcos0 
sin(2rrft) in phase vith the crystal motion which contributes to 
the inertial energy of the system, and a component Bsin 0cos 
(2rrft) 90 degrees behind the crystal motion which contributes 
to the damping of the crystal motion. Although damping is not 
included in the Rayleigh theory, its effect will be negligible for 
small deposit thicknesses, since sin 0 ~ 0 . At larger thick- 
nesses, however, its effect may increase to significant pro- 
portions, as noted in the Discussion. 


Droplet Deposit 


Figure lb shows the model of the motion of the oscillating 
liquid droplet. If surface tension forces are ignored, the motion 
is described by equation (11), but the boundary conditions are 
now: 

v = v sin 2rrft at x = 0, 


and 



0 at all points on the droplet surface 


where v is the tangential velocity on the droplet surface, and 
n is the coordinate normal to the surface. The solution to 
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equation (11) with the above boundary conditions will be much 
more complex than for the uniform film, because the variable 
droplet geometry requires the introduction of more descriptive 
parameters. A solution to this problem has not been located in 
the literature, and hence for the present purposes the velocity 
distribution in a droplet of height x will be assumed to be the 
same as that in a uniform film of tffckness x . In fact, the 
actual velocity at a gi'-en distance from the crystal surface 
she J be higher in the droplet than in the uniform film, because 
th ffect of area reduction with distance in the droplet will be 
to reduce the velocity gradient needed to exchange momentum 
with the outer levels of the droplet. The loss of QCM response 
calculated using the uniform film velocity profile will therefore 
be greater than the actual loss. 

DETERMINATION OF THE RESPONSE FACTOR 

It was noted earlier that the loss of sensitivity of the QCM 
to a liqu ; 1 deposit is due to three effects which act together but 
which ; e conceptually separable: viscous effects, droplet 
forma* on, and droplet distribution. In the interest of clarity of 
presentation, these three effects are considered one at a time in 
the following analysis. F is derived in turn for the uniform 
liquid film, individual liquid droplets and then distributed liquid 
droplets . 

Tn the theoretical analysis, F is shown to be primarily a 
function of the distance of the deposit from the QCM surface, 
while the primary data output from the QCM is miss per unit 
area. To connect the theory with experimental data, the concept 
of nominal deposit thickness, £ , is used. £ is defined as the 
deposit mass per unit area divided by the density. Further.it 
was shown in the viscous flow analysis above that distance from 
the QCM surface to points in the deposit can be non-dimensional- 
ized by multiplying them by k, equation (14). The theoretical 
expressions for F are thus developed as functions of k£ . The 
following paragraphs present the derivation of the kinetic energy 
of the various deposits, AT. Finally, in order to simplify the 
algebra, AT is derived per unit area, permitting the crystal area 
to be eliminated from the integrations. 

. riform Film Deposit 

For the uniform solid deposit the velocity throughout is 
equal to that of the crystal surface. The kinetic energy per unit 
crystal surface area, AT is given by 

.“V® 2 

AT = I (v sin 2rrft) dx 

SU J o ° (17) 

= ip x (v sin2nft) 2 
m o 
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By comparison with equation (10), the kinetic energy coefficient 
per unit area, Aa^, is given by 


Aa 


su 


0 * 


(18) 


The kinetic energy of the uniform liquid deposit, AT £ u wiU be 
a function of the component of the local velocity in phase with 
the 5 nertia forces, v B cos0sin2rrft. 


AT 


£u 



(v cos0 sin 2rrft) dx 
o 


(19) 



(B cost) )^ ’ dx (v 0 sin2TT f t)^ 


By comparison with equation (10) 


Aa 


£ u 



(B cos 0 dx 


( 20 ) 


From equations (9), (13) and (20) the response factor for a 
uniform film, F , is thus 


/■ m 2 

F =1 (B cos 0) dx/x 
u J m 

o 

= f (B cos 0 )“ d(x/x ) 

Jo m 


( 21 ) 


F has been computed from equations (12), (13) and (21) for 
various values of kx , and is plotted in Figure 3 as a function 

of k£ , noting that For a uniform film £ is equal to x . The 
, nu . h nu 1 m 

shape ot tr curve represents a transition from the very small 

response loss for small deposits (k£ <0.2) to the asymptotic 

response of a qua si -infinite film (k£ >1.6) which is almost 

*'tally confined to the region near the crystal and is thus inde- 
pendent of nominal thickness. Also plotted in Figure 3 is the 
product F • k£ , which is proportional to the actual QCM out- 
put. The rising of this QC.M output function to a maximum, 
followed by a steady decline is the most commonly encountered 
practical manifestation of QCM viscous effects. 

Single Droplets 

Whereas the geometry of a uniform film can be character- 
ized by a single dimension, several geometric parameters are 
necessary for the specification of condensed droplets. Since 
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the purpose of this paper is to demonstrate the credibility of the 

analysis rather than perform an analysis of droplet geometry 

some simplifying assumptions are made. The droplet is assumed 

to have the form of a spherical cap with contact angle 0 and 

height x . Relationships between area, volume and x c and 0 

j m , . . , . m c 

are derived in the Appendix. 

For a solid droplet of volume and base area the 

kinetic energy per unit bai 2 ar»a, AT^, is given by 


f T 


sd 




m 


A(x) (v sin2nft) dx/A 


= ip (V ,/A ) (v sin 2nft) 2 
aao 


( 22 ) 


Hence, from (A6), (A9) and (10), 

' 2 - 


Aa 


sd 


- ( 2 
"' r ' 


cos 0 (2+sin A ) 
c c 


(1-cosG ) sin e 
c c 


) 


(23) 


For a liquid droplet of the same dimensions, the kinetic energy 
per unit base area, aT^j* is (assuming that the velocity profile 
is given by equation (15)) 


AT 


Id 



A(x) • (v B cosf sin 2TTft) 2 dx/A 
o d 

A(x)/A.) (Bcos0) 2 dx* (v sin2nft) 2 
d o 


(24) 


From equations (10), and (24) 

x 


Aa 


Id 


= p/ 


(A(x)/A d ) (BcosO) dx 


(25) 


Hence from equations ( 9 ). (23) and (25) 


F d = 


3 sin 0 (l-cos0 ) 


(2-cos0 (2+sin 0 ) ) 
c c 


Z * 1 (A(xi 

J r\ 


)/A )(BcosG) d(x/x ) 
a m 

(26) 


In order to evaluate F^ for the purposes of illustration of the 
effect of droplet and/or comparison with the experimental data 
of the next section, it is necessary to know the value of 0 . 0 
depends on several parameters, such as temperature ancf c 
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surface condition, and data relevant to the present situation are 
not available. However, photomicrographs made by the author 
of similar deposits showed contact angles near 30°, and it can be 
shown by evaluating equations (A 5) and (A6) that the function 
Alxl/Aj varies little for 10° < 0_ < 50°. So that the derived 
equations may be evaluated a value of 0 equal to 30° has thus 
been assumed. F. has been computed frori^ equation (26) for 
various values oFkx using 0 equal to 30°, and equations (12), 
(13), (A5) and (A6). "Vhe nomuial deposit thickness of the drop- 
let, £ nt j relative to its own base area, A^ is given by 

'„d = (27) 

where V, and A are given by equations (A12) and (All). F^ 
is plotted versus ki , in Figure 4. While the relationship is 
qualitatively similar to that for the uniform film, the effect of 
configuring a liquid deposit as a droplet is shown to be a further 
decrease in the response of the QCM for a given nominal 
thickness. 

Distribution of Droplets 

A real liquid deposit will consist of a distribution of drop- 
lets of various sizes and spacings. The effect of droplet spacing 
will be to cause the deposit mass per unit area of QCM crystal - 
the quantity sensed by the QCM - to be lower than the mass per 
unit droplet base area. Assuming all droplets to be the same 
size, the nominal deposit thickness for the distributed droplets, 

£ , ,, is related to the nominal thickness of the droplet referred 
to its base area, by an area coverage factor, ar^: 


^ndd a A x ^nd 


(28) 


where 


a , 


area covered by droplets _ . N , 7Q , 
total QCM area d d 


The response factor, F , ,, is a function only of droplet geometry, 
not the distribution, ann^hus is given by equation (26) as a 
function of kx . However, the computed value of F^ is now 
plotted versus^! ,,, given as follows (from equations (27) and 
(28) ). ndd 

k£ aa ~ »a x 0. 512 kx (30) 

ndd a m 

Figure 4 shows a family of curves of F versus k£ , , for 
various values of a In practice has a maximum value of 
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tt/4 cos 30 or 0.907, corresponding to a close packed hexagonal 
distribution. Another significant value of or a * s 0. 785, for a 
close packed square distribution. The QCM response is seen 
to fall with decreasing area coverage. 

PRACTICAL APPLICATION OF THE THEORY 


The single fundamental output from a QCM with an added 
mass is change of frequency. A solid, evenly distributed 
deposit can be characterized by a single parameter-mass per 
unit area - and the QCM sensitivity is constant. Hence, a single 
frequency change data point can be translated into an added mass. 
However, unless the deposit mass is small enough for F to be 
near to unity, a liquid deposit requires several parameters to be 
determined, such as droplet shape and size and area distribution, 
and liquid viscosity. Hence several data points are required for 
determination of the deposit parameters, which precludes direct 
unequivocal measurement of the mass of a particular liquid 
deposit. The most apr- opriate use of the theory appears to be 
in the analysis of processes made up of a series of data points 
from which the various parameters can be deduced. Two 
situations of this type occur quite frequently in practice. In one 
situation a QCM with a deposit on its surface is warmed up with 
no loss or rearrangement of deposit. In this situation Ic varies, 
while the geometric parameters A^, V^, and a ^ remain 

constant. Thus the response factor F^ varies with k£ nc id at 
constant a ^ , and the QCM response during warm-up should 
follow one of the constant a A characteristics of Figure 4. In 
the other situation the QCM is held at constant temperature while 
the accumulation or loss of deposit is measured, in which case k 
is constant, while the geometric parameters vary. It is suggested 
that during condensation a liquid droplet deposit builds up by 
growth of the individual droplets at constant distribution density, 
N^, until adjacent droplets touch and coalesce, thereby reducing 
and increasing the spacing between droplets. Growth resumes 
at constant until coalescing occurs again. In practice, the 
droplets wi\l have a distribution of sizes and spacings, so a suc- 
cession of discrete steps in Nj are unlikely. A more likely 
process is growth at constant N,} from an initial condition of 
relatively widely and randomly spaced nucleation centers, until 
the droplets begin to coalesce, whereafter adjusts continuously 
downwards while the deposit maintains a more or less close 
packed arrangement. Characteristic paths can be derived for 
this postulated process as follows. Equations (28) and (29) can 
be combined to give 


k£ , , -• A • N ’k£ 
ndd d d nd 


(31) 


For 0 c equal to 30°, noting (27), (31) and (A6) 
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(32) 


“ndd = 167 <V k ‘ > 


(k *nd> 


The relationship between F and k£ for isothermal growth 
at constant N /k^ is found aa by determining F^d for a given 
k£ , from equations (26) and (27), andkf ,, from equation (32) 
fofconstant values of Nd/k^. These conJVant Nd/k^ charac- 
teristics are shown in Figure 4. According to the suggested 
mechanism growth will occur along one or more lines of con- 
stant Nd/k^ until a close packed situation occurs, at which 
point Nd will decrease and growth will continue along another 
line, constant Nd/k^. 


EXPERIMENTAL INVESTIGATIONS 


The validity and usefulness of the theory presented above 
was examined by a series of experiments using Dow Corning 
DC 704 silicone oil on a Celesco Model 700 quartz crystal 
microbalance. Noting the remarks in the previous section, two 
types of experiments were conducted. In the first type the 
response of constant temperature QCM to a steadily increasing 
deposit was measured. In the second type a deposit was formed 
on the QCM and then its response was measured as a function 
of temperature. 


The apparatus is shown schematically in Figure 5. The 
QCM can be maintained at any desired temperature above 77°K 
by balancing the electric power into a resistance heater wound, 
on the- holder against heat lost through a thermal link between 
holder and shroud. Also mounted within the shroud but shielded 
from the QCM by a cooled partition is a pot containing DC 704. 
The pot temperature can be maintained at any temperature above 
77°K by balancing electrical heater pov/er against a thermal i. 

The pot has an orifice whose axis lies along the line between pot 
and QCM permitting a flux of DC 704 molecules to be directed 
at the QCM sensing crystal through a hole in the partition, which 
can be closed by a remotely controlled shutter. The whole 
apparatus is mounted in a glass bell jar, evacuatable to pressures 
below about 2x1 0"? torr. 

Viscosity of DC 704 

In order to determine k for an experimental situation, 
the relationship between the viscosity of the subject liquid must 
be known as a function of temperature. Figure 6 shows kine- 
matic viscosity versus temperature obtained from the Dow Corn- 
ing Company. DC 704 is a phenyl-methyl siloxane, a high 
molecular weight silicone oil related to the general family of 
polydirnethyl siloxanes, which have been shown to exhibit non- 
Newtonian flow behavior at high shear rates. (5) The shear 
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rates experienced by a liquid deposit on a oscillating QCM can be 
quite high, and could conceivably reach the order of magnitude 
required for non-Newtonian behavior to occur. The possibility 
of this situation being reached in the present experiments is 
covered in the Discussion. 


QCM Response Data 

In the first series of tests the QCM was held at constant 
temperature while its surface was exposed to a constant flux of 
DC 704 molecules. The test was begun with the QCM surface 
clean. The shutter was then opened and the QCM output re- 
corded as a function of time. Data interpretation was based on 
the assumption that after the initial adsorption and nudeation 
processes were completed, the condensation rate could be as- 
sumed to be constant. This constant actual condensation rate 
was determined from the initial indicated condensation rate, 
since the response factor for small nominal thicknesses is es- 
sentially unity. The actual mass per unit area on the QCM at 
later times in the deposition was then calculated from the assumed 
constant condensation rate and the time since the beginning of 
exposure. The accuracy of this linear extrapolation was con- 
firmed by interrupting the flux and cooling the QCM at the con- 
clusion of the test in order to regain full response of the deposit 
by solidification. The agreement between this final measurement 
and the magnitude predicted by the linear extrapolation was 
surprisingly good -- always better than five percent, and usually 
no worse than one percent. The response factor was calculated 
as a function of deposit thickness by dividing the indicated mass 
per unit area by the calculated mass per unit area for the same 
instant of time. The nominal thickness was found from the cal- 
culated mass per unit area and the density, while k was found 
for the test temperature from equation (14) and Figure 6. Data 
from four experiments of this type, performed for QCM temper- 
atures of 27°C, 37°C, and 51°C, are plotted in Figures 7 and 8. 
The experiment at 37°C was repeat* r or two net condensation 
rates -- 0. 61 x 10"® and 1. 22 x 10"® gms/cm^/ sec. The data 
for 27°C and 51°C were obtained for net condensation raves of 
1. 53 x 10"® and 1. 72 x 10 _ 8 gms/cm^/ sec, respectively. The 
data for 27°C and 37°C appear to follow the scenario proposed 
in the previous section. For the early stages of growth (k£ n < 
0.25) the data appear to follow roughly, two Nj/k^ character- 
istics, with a transition near k£ ~0. 1. For 0.25 <k£ <0.5 
the process occurs along a more or less close-packed character- 
istic. The departure of the data fr 'm the theoretical close packed 
characteristics for k£ n > 5 is prouably due to the increasing 
significance of effects neglected in the present analysis, and is 
covered in the Discussira. The number of droplets per unit area 
in these deposits can be estimated by fitting theoretical constant 
N d /k2 curves to the data for 0.05 < kX n <0.25, using equation 
(31) or Figure 4. 




At 37°C, for which k is 11687 cm' , the data can be 
fitted by N<j/k^ equal to about 0.045, making of the order 
of 6x10® cm"2. At 27°C, for which k is 9474 cm"*, N^/k^ 
equal to 0.04 is a reasonable fit, making about 4x10 6 cm"2. 
Use of these values for at 51°C, where k is equal to 14980 
cmkndicates that corresponding values of N^/k^ at 51°C should 
be in the range of about 0.03 to 0.02. Constant Nj/k^ curves 
for these values have been drawn on Figure 8, showing that 
51°C data do fall in this range, even though the shape of the 
data is inconclusive. 

These values of Njj can be compared with the photo- 
micrographs of Shapiro and Hanyok'®*, which show DC 704 
deposits on various surfaces. These photomicrographs show 
number densities in the range of 5x10* to 1x10° cm"‘ for 
deposits on a mirror surface, which, because of its smooth- 
ness would be expected to have somewhat fewer nucleation sites 
than the QCM surface. 

In the second series of tests deposits of different 
magnitudes were condensed on the QCM and their response 
measured as a function of temperature. The true mass per 
unit area of each deposit was found by cooling the QCM below 
the solidification temperature, which was near -35°C. The 
response at other temperatures was determined by heating or 
cooling the QCM through a series of equilibrium temperatures. 
This technique is limited to temperatures below the region in 
which re-evaporation rate of the DC 704 became significant, 
which was about 20 C. The response factor at a given tempera- 
ture was calculated from the mass per unit area indicated by the 
QCM at that temperature divided by the true mass per unit 
area. The nominal thickness was found from the actual mass 
per unit area divided by the density, k was variable in these 
tests and was calculated from equation (14) and Figure 6. The 
data for several experiments with different nominal thicknesses 
are shown in Figure 9. The maximum value of k£ reached in 
each test was determined by either the attainment n of a tempera- 
ture at which the re -evaporation rate became significant (about 
20°C) or an abrupt loss of output occurred, the possible cause 


of which is suggested in the Discussion. It was noted in the 
previous section that the deposit geometry should remain con- 
stant during these tests, in which case the data should follow 
the constant ota characteristics of Figure 4. The data do indeed 
follow the same general pattern as the theoretical curves, 
although they appear to show a slightly faster decline at high 
values of k£ n . However, the data for the 5.12x10“^ cm deposit 
shows the same inflection point near k£ n of 0.4 as predicted by 
theory, this being the only set of data in which the inflection 
point was observed. 

The slightly different shape of the data curves from the 
constant a A curves of Figure 4 is thought to be due to tempera- 
ture dependent effects ignored in the theory, for example, drop- 
let shape or constant angle. It was first thought that the 
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discrepancy might have depended on the temperature at which 
the deposit was formed. A deposit formed on a cold QCM may 
have a continuous film morphology, which breaks down into 
droplets as heating occurs. To investigate this possibility 
several experiments were made with a 6.54 x 10"^ cm film 
thickness. In one experiment the deposit was formed cold 
(-45°C) and the data were acquired during heating. In a second 
experiment an equivalent deposit was formed near 20°C and 
data were obtained during cooling, then by reheating and re- 
cooling. The data from these two series of tests shown in 
Figure 9 are virtually indistinguishable. 

It wao concluded from these data that a transition from 
continuous film to droplet geometry, or any other irreversible 
geometric change does not occur in the temperature range 
studied, even if the deposit is formed in the solid phase. This 
means that the surface mobility of DC 704 at -45°C is sufficient 
to produce a morphology similar to that obtained for liquid 
phase deposition. 

The deposit nominal thicknesses used in these tests was 
quite high, so that sufficiently large values of k£ could be 
reached. It is possible tint the thicknesses were high enough 
to cause the deposit morphology to begin to differ from the 
distributed, even sized droplet pattern assumed with apparent 
success for smaller nominal deposit thicknesses. This pos- 
sibility was assessed using the following relationship, which 
can be derived from equations (27), (29), (30), and (A6). 


y =5. 
a 


5- i 


2/3 


ndd 


N 


1/3 


(33) 


,-5 


The data for JJ ^ equal to 1.75 x 10 cm can be fitted 
quite well by the theoretical curve for a ^ equal to 0. 6, 

Figure 4. Insertion of these numbers in equation (33) gives 
equal to 4.2x10°- 2 cm which is in good agreement with the 
earlier data, and tends to confirm the validity of the model for 
this nominal thickness. The data for 5.12 x 10 - 5 crn for 
k^ndd can fitted by a A equal to about 0.8, while for 

k£ n dd >0.25 a better fit is obtained for a \ equal to 0. 70. 

These y A values correspond to Nj values of 1.2x10^ -2 cm and 
8.0 x 105 -2 cm , respectively, and also suggest that at this 
thickness the deposit is approaching a close packed situation. 
The data for the 6. 54 x 1 0-b C m deposit at kjt ^ . < 0. 3 can only 
be fitted by y ^ equal to unity, which is impossible for a droplet 
deposit, for which yj^ cannot exceed 0.91. This seems to 
indicate that a transition to a more film -like morphology is 
occurring at these higher nominal thicknesses, which is to be 
expected. 
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DISCUSSION 


The intent of this paper was to demonstrate that the 
variable response of the QCM to liqui deposits can be ex- 
plained almost totally by viscous effects coupled with a droplet 
type deposit morphology. It is apparent irom he agreement 
obtained between the theory and the experirm nts that this has 
been satisfactorily achieved. Good qualitative and quantitative 
agreement was obtained up to non-dimensional nominal thick- 
nesses of about 0. 5. Above this value a very rapid fall off in 
response was experienced in all tests where this k£ value was 
reached. Detail analysis of the experimental data in this region 
revealed a discontinuity in the second differential of F versus 
k£ R , indicating a change of controlling phenomenon. It is sug- 
gested that this could be due to the onset of non-Newtonian 
effects in the silicone oil, or a change in the modal response of 
the oscillating crystal under the higher damping forces which 
will arise with heavier deposits ^ .The oscillation amplitude of 
the crystal is of the order nf 50A,''' so that for a deposit of 
nominal thickness of 1000 A on a 10MHz crystal the shear rate 
is of tlie order of (1 0?x2x50)/ 1 000 or 10^ sec" . This is well 
within the range in which silicone fluids can show a sharp loss 
of apparent viscosity. However, a detailed analysis of the 
velocity profile in the deposit, plus viscosity data for DC 704 
at high shear rates would be needed to investigate this point 
further. Assessment of the possibility of modal changes is 
beyond the scope of this work. 

Although the quantitative agreement between tne theory 
and the data appears to be quite good, it should be re- 
emphasized tnat several major assumptions were made. 

These were: 

(1) Use of the uniform film velocity profile for the 
liquid droplet. 

(2) Assumption of spherical -cap droplet shape and 
uniform droplet size distribution. 

(3) Neglection of damping effects in the Rayleigh 
frequency perturbation analysis. 

(4) Exclusion of all other possible surface forces 
except viscous drag from the analysis. 

(5) To permit reduction of the data the viscosity data 
obtained from Dow Corning were arbitrarily 
extrapolated from about -12°C to about -35°C. 

It is logical to suppose that refinements of the analysis 
in these areas would produce better agreement with the 
experiment. 

It was noted in the text that it is not possible to use the 
theory to interpret isolated data points because of the many 
parameters involved. Instead the utility ol the theory will 
probably be to generate insight into condensation and revaoora- 
tion processes. Interpretation of the data presented in this 
paper has ger. erated credible values for area coverage and 
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number density for the droplet and has suggested that the 
spherical cap shape assumption may be acceptable. Further 
work in which some of the neglected effects are included in the 
analysis, coupled with more careful, systematic experimental 
measurements may well prove to be a fruitful source of data 
on contamination deposit morphology. 

NOMENCLATURE 


d 

a, Aa 
B 

c, Ac 
F 
f 
k 

(A rn / A ) 

N d 

r 

S 

T 

u 

V 

V d 

v 

Vr, 


S A 

Ih 

o 


Droplet base area 

Rayleigh kinetic energy coefficient 
Oscillation velocity magnitude 
Rayleigh potential energy coefficients 
QCM response factor 
QCM frequency 
(2nfp / 2^)z 

Deposit nominal thickness 
Deposit mass per unit area 
Number of droplets per unit area 
Solid radius of spherical cap droplet 
QCM sensitivity 
Rayleigh total kinetic energy 
Rayleigh displacement coordinate 
Rayleigh total potential energy 
Droplet volume 
Local velocity in deposit 
Maximum crystal surface velocity 
Distance between QCM surface and a point 
in the deposit 

Distance between QCM surface and the most 
distant point in the deposit 
Fraction of QCM surface covered by droplets 
Phase angle of motion in the deposit 
Droplet contact angle 
Kinematic viscosity of the deposit 
Deposit density 


Subscripts : 

i droplet 

dd distributed droplet 

l liquid 

n nominal 

s solid 

u uniform 
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APPENDIX 

The principal geometric parameters of the spherical cap 
droplet are indicated in Figure 10. Other parameters may be 
derived as follows: 

Droplet height, x: 

x = r(cos9 - cose ) (Al) 

c 

x = r(l - cose ) (A2) 

m c 

Hence x = x (cos© - cose )/(l - cose ) (A3) 

m c c 


Droplet cross-sectional area, A(x): 

2 

A(x) = TT(rsine) 

2 2 

Obtaining r from A2 and sin 6 from A3 


(A4) 


A(x) = 


TTX 


(i - cos e ) 
c 


1 - ( (x/x ) (1-COS0 ) +COS0 ) 
m c c 


(A5) 


259 


Base area, A , (= A(x) at x/x = 0): 
d m 


nx ^ . 2„ 

A = m sm e c 

“ 2 
(1 - cosG^F 


(A6) 


Volume, V,: 
d 


r m 

A(x) • dx 


(A7) 


Substituting A4 in A7 


V d = 


nr 


2 - cosO c (2 + sin e c )j (AS) 


Substituting A 2 in A8 


V , = 


2-cosG (2 + sin^p ) 
c c_ 

(1 - cose 

c 


(A9) 


For G c equal to 30° the following relationships follow from 
A5, A6 and A9. 


A(x) = 175 x 


l-(0. 134(x/x )+ 0,886) 

m 


Aj = 43. 8x 
d m 


(A10) 

(All) 


V , = 22.4 x J 
d m 


(A12) 
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SURFACE 


(a) OSCILLATING UNIFORM I IQUID FILM DEPOSIT 



(b) OSCILLATING LIQUID DROPLET DEPOSIT 


Figure 1. Schematic of Uniform and Droplet 
Liquid Deposits. 
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POSITION IN OSCILLATING FILM, 

Vx,„ 

Figure 2. The Dependence of Amplitude, B , 
and Phase Angle 0, on Pos : t>on in the Liquid 
Film, x/x m and Non-Dimensional Thickness, 
kx m , for an Oscillating Liquid Film. 



RESPONSE FACTOR 



NOMINAL THICKNESS. kf n 


Figure 3. Theoretical Response Factor and QCM 
Output Versus Deposit Thickness for Uniform 
Film Liquid Deposit. 
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Fiaure 4. Theoretical Relationship Between 
Response Factor and Nominal Thickness for 
Liquid Droplet Deposits = 30°). 
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RESPONSE FACTOR, 



NON-DIMENSIONAL THICKNESS, k I 

n 

Figure 7. Experimental Data for Response 
Factor as a Function of Deposit Magnitude 
at 37°C. 



NON-DIMENSIONAL THICKNESS, kf 

n 

Figure 8. Experimental Data for Response 
Factor as a Function of Deposit Magnitude 
at 27°C and 5T°C. 
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NON-DIMENSIONAL THICKNESS, k » n 
Figure 9. Experimental Data far Response 
Factor as a Function of kf n for Three Values 
of Nominal Thickness. 
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Figure 10. Major Dimensions of a Spherical 
Cap Droplet on a Surfoce. 


